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ABSTRACT 


'"in  implicit  numerical  model  for  two-dimensional  hydrody¬ 
namic  flow  in  coastal  seas  by  Leendertse  (1967),  as  modified 
by  Hart  (1976)  ,  was  applied  to  Monterey  Bay.  The  model  was 
tested  against  available  water-level  and  current  observa¬ 
tions.  The  responses  of  Monterey  Bay  to  tidal  forcing  and 
steady-state  winds  were  simulated.  Under  tidal  forcing  it 
was  found  to  provide  reasonable  estimates  of  sea-surface 
elevations,  currents  were  not  well  predicted,  indicating 
that  other  mechanisms  such  as  wind,  density  stratification, 
and  oceanic  currents  generally  dominate  the  forcing  of  the 
circulation  in  Monterey  Bay.  The  andel  in  its  present  form 
was  found  to  be  potentially  suitable  for  providing  real-time 
tid9  correctors  during  a  hydrographic  survey,  achieving  an 
BBS  error  of  H.5  cm  in  predicting  sea-surface  elevations 
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This  study  grew  out  of  a  desire  to  extend  tidal  data 
observed  at  a  fev  locations  to  the  antire  area  of  a  hydro- 
graphic  survey,  with  sea-surface  alavations  modeled  c  ar 
the  whole  field*  survey  deaths  aay  be  corrected  autoa  L- 
cally  by  subtracting  realistic  values  of  the  surface*: 
variation  froa  datus  at  any  point*  at  any  tiae. 

Tidal  zoning  to  obtain  sea-surface  elevations  is*  - 
present*  a  subjective  affair  requiring  numerous  water-level 
stations  to  indicate  the  progress  of  a  tidal  wave  into  an 
inlet.  The  pattern  of  propagation  at  points  distant  from 
the  observing  stations  is  inferred  only  qualitatively  froa 
bathyeetry.  Correctors  are  determined  by  defining  zones 
graphically  and  computing  appropriate  phase  and  amplitude 
adjustaents  for  each  zone  to  apply  to  tidal  values  observed 
at  the  reference  stations.  This  process  does  not  provide  a 
continuum  of  correctors.  It  requires  subjective  judgment 
and  considerable  experience  to  achieve  adequate  results. 

The  analysis  is  typically  performed  well  after  the  survey* 
when  it  is  too  late  to  use  depth  data  corrected  for  observed 
tides  to  provide  cross-checks  on  tha  positional  accuracy  of 
rhe  data.  2rrors  that  might  have  been  detected  and 
corrected  during  the  survey  nay  pass  unnoticed  until  an 
expensive  return  to  the  survey  area  or  downward  classifica¬ 
tion  of  the  survey  is  necessary. 

Ose  in  the  field  of  a  two-dimensional*  numerical  model 
for  circulation  and  sea-surface  elavation  would  alleviate 
these  difficult ies.  Such  a  model  must  be  simple  and  flex¬ 
ible  if  it  is  to  be  applied  in  real  or  near-real  time  on 


field-type  microprocessors.  In  coastal  seas  and  inlets,  it 
should  be  able  to  provide  sufficiently  accurate  surface 
elevations. 

To  test  this  concept,  the  two-dimensional,  hydrodynamic 
model  of  Leendertse  (1967)  ,  as  modified  by  Hart  (1976)  and 
during  this  study,  was  applied  to  Monterey  Bay.  The  model 
had  previously  bean  used  with  some  success  in  shallow 
coastal  seas  and  estuaries,  but  not  in  an  area  with  such  a 
dominant  bathymetric  feature  as  the  Monterey  canyon. 

Although  no  attempt  was  made  to  compare  this  model  to  any  of 
the  broad  spectrum  of  models  in  use  throughout  the  oceano¬ 
graphic  and  coastal-engineering  communities  [Tracor,  Inc., 
1971],  its  relative  simplicity,  ease  of  implementation, 
flexibility,  and  accurate  output  are  all  important  to  its 
potential  usefulness  as  a  tidal  zoning  system  during  field 
operations. 

An  additional  purpose  of  this  study  is  to  incorporate 
large-scale  non-tidal  forces  into  the  model  to  explore  their 
effects  on  the  circulation  and  sea-surface  elevation  of  a 
coastal  body  of  water.  Various  investigators  have 
suggested,  in  fact,  that  tidal  forces  are  overridden  in 
their  effect  on  the  circulation  of  Monterey  Bay  by  the 
influence  of  offshore  currents  and  atmospheric  conditions. 
The  potential  significance  of  such  factors,  and  the  ability 
to  incorporate  real-time  observations  of  them  into  the 
model,  are  also  important  to  the  model's  usefulness  as  a 
tidal  zoning  system. 

B.  HISTOBT  OP  THE  HODEL 

The  original  version  of  the  numerical  model  used  during 
this  study  was  described  by  Leendertse  (1967)  .  His 
"multioperational"  finite-difference  scheme,  using  both 
implicit  and  explicit  techniques  to  solve  the  equations  of 
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fluid  «otion,  provided  advantages  in  computational  stability 
and  efficiency  over  the  explicit  models  then  current  [Hart, 
1976  ].  In  particular,  the  aodel  retains  stable  regardless 
of  the  time  step  used;  in  a  relatively  deep  embayment  such 
as  Honterey  Bay  the  investigator  is  not  restricted  to  time 
steps  of  the  order  of  seconds  as  is  common  with  explicit 
models  such  as  that  used  by  Lazanoff  (1971)  in  his  study  of 
the  bay. 

The  model  has  been  widely  applied,  both  in  small  harbors 
[Leendertse,  1967;  Leendertse  and  Lu,  1975;  Chiang  and  Lee, 
1982]  and  in  coastal  seas  [Leendertse,  1967;  Hart,  1976; 
Spaulding  and  Beauchamp,  19  83].  in  past  applications, 
length  scales  ranged  to  290  km  and  depths  ranged  to  100  m. 
This  study  extends  the  aodel  to  a  much  deeper  area  with 
pronounced  vortical  relief. 

C.  CIRCULATION  STUDIES  OF  HOITBREY  BAY 

Honterey  Bay  is  a  relatively  large  (16  by  42  km),  nearly 
symmetric  embayment  on  the  central  coast  of  California.  Its 
most  notable  bathymetric  feature  is  the  Monterey  canyon, 
which  curves  into  the  bay  from  the  southwest,  severing  the 
continental  shelf.  Within  the  bay  proper,  depths  rise  from 
750  m  at  the  seaward  end  of  the  canyon  to  an  average  55  m  on 
the  shelf. 

The  bay  presents  to  the  Pacific  Ocean  one  uninterrupted 
open  boundary,  some  36  km  long.  Consequently,  oceanic  tides 
and  currents  and  offshore  atmospheric  effects  are  primary 
driving  forces  of  circulation  within  the  bay  itself.  Local 
winds  and  seasonal  river  runoff  may  have  some  effects,  espe¬ 
cially  in  the  shallower  portions  of  the  bay  to  the  north  and 
south.  The  relative  importance  of  these  various  influences 
is  not  well-understood. 
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Figure  1.1  Hap  of  the  Study  Area,  Showing  Hod el  Grids 


Figure  1.2  Bathymetry  of  Honterey  Bay  Viewed  from  southwest 


Three  oceanographic  "seasons"  for  this  portion  of  the 
coast  were  first  described  by  sltogsberg  (1936).  From 
November  through  Pebruary  the  Davidson  Current  flows  north¬ 
ward  along  the  coast  in  conjunction  with  southerly  or  weak 
northward  winds  and  the  onshore  transport  of  surface  water. 
Prom  Harch  through  August  a  period  of  upwelling  is  accompa¬ 
nied  by  the  southward  flowing  California  current,  strong 
northwest  winds,  and  offshore  transport  of  surface  water. 

In  September  and  October  the  oceanic  period  is  marked  by 
relative  calm  and  an  increase  in  surface  temperatures. 

All  available  data  concerning  circulation  within  the  bay 
proper  were  summarized  and  analyzed  in  1973  as  part  of  a 
major  oceanographic  study  [Scott,  1973].  Normal  circulation 


was  found  to  be  northward  through  the  bay,  with  a  small  gyre 
forming  in  the  southern  bight,  A  more  recent  analysis  by 
Broenkow  and  Saethie  (1978)  also  concluded  that  flow  is 
generally  northward,  with  water  entaring  primarily  along  the 
axis  of  the  canyon  and  having  a  residence  time  of  2  to  14 
days  during  upwelling  periods.  They  also  suggest  that  a 
volume  of  10»  m3  pumped  into  and  out  of  the  bay  by  internal 
tidal  mixing  may  be  responsible  for  the  frequent  presence  of 
cool,  nutrient-rich  waters  at  the  head  of  the  Monterey 
Canyon  near  Moss  Landing. 

Most  circulation  studies  of  the  bay  proper  have  relied 
primarily  on  temperature  and  salinity  data  collected  at 
oceanographic  stations.  Although  drift  cards  and  similar 
devices  have  been  used  to  map  surface  currents,  long-term 
current-meter  observations  have  not  been  available  until  the 
last  ten  years.  During  pre design  studies  for  sewer 
outfalls,  current  meters  were  deployed  for  periods  of  a  year 
or  more  near  Santa  Cruz  [Brown  and  Caldwell,  Inc.,  1978  ], 
the  Pajaro  River  [Environmental  Research  Consultants,  Inc., 
1976  ],  and  the  Salinas  River  [Engineering-Science,  Inc., 
1977],  In  general,  these  studies  suggest  that  th§  net  flow 
of  water  is  northward  along  the  coast,  with  some  dependence 
on  the  local  diurnal  wind. 

Tidal  forcing  has  not  been  examined  closely  in  any  of 
the  aforementioned  studies.  However,  an  explicit  numerical 
modal  was  employed  by  Lazanoff  (1971)  to  study  the  tidal 
circulation  within  the  bay.  Although  field  observations 
indicated  that  the  primary  driving  force  for  circulation 
derived  from  oceanic  currents,  the  tide-  and  wind-driven 
model  did  predict  correct  sea-surface  alevations  and  current 
phases  and  directions  for  the  short  time  periods  over  which 
it  could  be  run.  Current  magnitudes  appeared  too  large  near 
coastal  boundaries. 
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More  recently  Br  etschneider  (1980)  analyzed  the  effects 
of  various  oceanographic  conditions  on  the  sea-level  changes 
observed  at  Honterey.  Variations  in  geostrophic  current 
floe,  ateospheric  pressure,  sea-surface  tenperature,  and 
eeridional  wind  stress  were  shown  to  correspond  to  observed 
variations  in  sea-surface  elevation  at  Monterey. 
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ii.  zsss&imss  si  in  laasneu  asfisi 


A.  HTOBODTIAHIC  THEORY 

The  numerical  aodel  of  Leendertsa  (1967)  ralies  on  the 
basic  equations  that  describe  conservation  of  momentum  and 
aa3s  for  incoapressible  fluid  notion.  In  the  Cartesian 
coordinate  system  of  the  aodel,  vith  x-  and  y-axes  eabeddad 
in  a  horizontal  f-plane  tangent  at  the  origin  to  the  undis¬ 
turbed  sea  surface  (the  datum)  and  with  the  z-axis  oriented 
upward,  these  well-known  aquations  are: 


5u  6u  5u  j5u  1  5p 

St  +  u5x  +  vS7  +  w5I  =-pS7+Fx 


5v  <5v  <Sv  Sv  1  dp  .  „ 
St  +  US7+VS7+WSI  “  "p  3*  Fs 


6w  6w  5w  1  dp  ,  _ 

+  +  *  "!*  '  —  J*  *  F! 


du  5v  dw 

S?  +  S7  S? 


(2.1) 


(2.2) 


(2.3) 


(2.4) 


Tha  variables  u,  v,  and  w  are  components  of  velocity 
parallel  to  the  x-,  y-,  and  z-axes  raspectively,  p  is  pres¬ 
sure,  and  p  i3  density.  Ihe  applied  forces  per  unit  mass 
(P  )  represent  effects  of  the  Earth’s  rotation,  the  Earth's 
gravitation,  viscous  and  turbulent  stresses  in  the  fluid, 
and  astronomical  tides. 
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These  equations  are  siaplified  by  malting  assumptions 
appropriate  to  the  examination  of  loag-period  forcing  in  a 
two-dimensional,  shallow  field.  Detailed  derivations  may  be 
found  in  Leendertse  (1967)  and  Hart(1  976).  &  more  general 

development  of  shallow  water  equations  may  be  found  in 
Csanady  (1982).  The  necessary  assumptions  are  summarized  in 
the  following  paragraphs. 

First,  because  a  coastal  sea  or  estuary  is  generally 
shallow  relative  to  the  horizontal  scale  of  motion,  the 
vertical  velocity  is  assumed  small  relative  to  the  hori¬ 
zontal  velocities.  Therefore,  both  convective-inertia  terms 
and  rotational  effects  that  involve  the  vertical  velocity  in 
equations  2.1-2.  3  may  be  neglected. 

Second,  the  equations  2.1-2. 4  are  averaged  to  model 
fluid  motions  with  periods  greater  than  those  of  short- 
period  turbulent  motions. 

Third,  the  hydrostatic  approximation  is  made  by  analysis 
of  equation  2.3.  The  vertical  component  of  the  rotational 
effect  may  be  considered  negligible  (of  the  order  10-* 
cm/s2)  and  so  may  vertical  stress-gradient  effects  (10- 2 
cm/s2  according  to  Csanady,  1982)  .  Furthermore,  since  mean 
vertical  velocties  are  unlikely  to  be  greater  than  10  cm/s, 
over  sufficiently  long  time  periods  ( >1 03  s  or  15  minutes) 
the  total  vertical  acceleration  will  be  of  similarly  small 
order.  Neglecting  for  the  moment  tidal  effects,  an  expres¬ 
sion  for  pressure  may  be  obtained  by  integrating  the 
remaining  terms  over  depth: 

n 

P  a  P  +  9  p<5z  (2.5) 

2L 

Z 

In  this  expression,  n  is  the  sea-surface  elevation  and  pa  is 
the  atmospheric  pressure  at  the  sea  surface,  both  functions 
of  x  and  y.  The  gravitational  acceleration,  g,  is  assumed 
constant  and  equal  to  its  mean  value  at  the  undisturbed  sea 
surface. 
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Equation  2.5  permits  computation  of  pressure  gradients, 
5p/6x  and  <5p/<$y,  froa  horizontal  gradients  in  sea-surface 
elevation  and  density.  Gradients  in  atmospheric  pressure 
aay  be  neglected  since  their  effect  is  small  relative  to  the 
turbulent  stress  induced  at  the  surface  by  the  wind.  In  the 
numerical  modal,  tidal  forcing  is  accomplished  by  generating 
gradients  of  sea-surface  elevation  rather  than  by  attempting 
tc  simulate  directly  the  astronomical  forces  that  cause  the 
tides. 

Fourth,  the  Boussinesq  approximation  is  made,  in  which 
the  influence  of  vertical  density  variations  is  assumed  to 
be  negligible.  For  this  to  be  true  the  area  to  be  modeled 
must  be  vertically  well-mixed,  an  assumption  that  is  not 
generally  applicable.  Although  the  version  of  Leendertsa's 
model  used  in  this  study  requires  this  assumption,  some 
compensation  for  density  stratification  might  be  made  by 
modifying  the  model  to  integrate  estimated  values  for  the 
vertical  density  variation  over  the  depth  at  each  point  to 
obtain  the  additional  density- induced  sea-surface  elevation 
[Csanady,  1982]. 

Fifth,  the  mean  viscous -shear  stresses  of  the  fluid  are 
assumed  negligible,  leaving  only  the  turbulent  stresses 
(Beynolds  stresses)  at  boundaries  within  and  external  to  the 
fluid  to  be  formulated.  Of  these,  3harp  density  gradients 
within  the  fluid  are  neglected  as  a  source  of  stress. 

Closed  lateral  boundaries  are  considered  by  applying  the 
coastal  boundary  condition  that  velocity  into  the  boundary 
is  zero.  Two  other  boundaries  are  considered,  the  sea 
surface  and  the  bottom,  and  algorithms  for  modeling  stresses 
on  these  boundaries  must  prepared. 

Sixth,  since  interior  stresses  resulting  from  sharp 
density  boundaries  are  assumed  negligible,  equations  2.  1  and 
2.2  aay  be  integrated  vertically  to  provide  implicit  expres¬ 
sions  for  horizontal  velocities  averaged  over  depth. 
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Applying  the  kinematic  boundary  condition  at  the  free 
surface  and  at  the  bottos  (assumed  impermeable),  equation 
2.4  say  also  be  integrated  vertically  using  the  Leibnitz 
rule  of  integration.  Three  equations  implicit  in  three 
unknowns  are  then  available  for  the  vertically  integrated, 
horizontal  velocities  0  and  V  and  the  free,  sea-surface 
elevation  n  : 

H  *  “ET  *  VJ7  -  -9  ♦  tv  *  (FMx  -  FBlt>/p(h+n>  <2-6' 

Hf  +  *57  *  '5  §■  ‘  *D  +  <FWy  '  FByl/p(h+,')  (2.7) 


5n  .  6(h+n)u  .  <S(h+n)v  „ 

_+ — — sT~-°  11 

In  these  equations  f  is  the  Coriolis  parameter,  h  is  the 
depth,  pwiis  the  surface  friction  stress  due  to  wind,  and  ? 
is  the  bottom  friction  stress.  0  and  7  are  mean  horizontal 
velocities  over  the  water  column.  This  simplification 
results  in  a  two-dimensional  model  with  which  patterns  of 
circulation  and  sea-surface  elevation  may  be  quantitatively 
determined. 

Finally,  3ince  bottom  stress  depends  on  the  fluid 
velocity,  a  well-known  quadratic  model  is  assumed  so  that 
the  stress  term  may  be  incorporated  directly  into 
Leendertse*  s  computational  scheme.  The  formulation  for 
bottom  stress  is: 

fbx  58  pgo(u2+v2)Vc2 


(2.9) 


f  »  -  pgV(02+V2)>5/C2  (2.10) 

By 

The  eepirical  Chezy  coefficient,  C,  may  be  coeputed  in  any 
of  various  ways  and  aust  be  specified  for  the  area  to  be 
modeled. 

B.  STBDCTOBE  AHD  COBPOBEHTS 

k  derivation  of  the  numerical  model,  a  discussion  of  its 
computational  stability,  and  a  program  listing  are  given  by 
Leendertse  (1967).  Key  features  of  the  computational  scheme 
and  the  computer  program  used  during  this  study  are 
presented  here. 

1.  computational  SSfr&»g 

In  the  numerical  model,  equations  2.6-2. 8  are 
approximated  with  a  finite- differencing  scheme  extending 
over  two  time  levels,  each  a  half  time  step.  At  the  first 
half  time  step,  t+1/2,  the  velocity  ff(tt1/2)  and  the  sea- 
surface  elevation  n(t+1/2)  are  computed  implicitly  and  the 
velocity  V(t+1/2)  is  computed  explicitly.  At  the  second 
half  time  step,  ttl,  the  velocity  V(t+1)  and  sea-surface 
elevation  n(t*1)  are  computed  implicitly  and  0(t+1)  is 
computed  explicitly. 

Computations  are  spatially  controlled  by  a  uniform 
grid  of  squares  laid  over  the  f-plaae  (Figure  2.1).  Depths 
relative  to  the  undisturbed  sea  surface,  here  taken  to  be 
mean  lower  low  water  (HLLW)  ,  must  be  supplied  for  the 
corners  of  each  square,  values  of  horizontal  velocity  are 
computed  at  the  centers  of  the  sides  of  each  square,  and 
values  of  sea-surface  elevation  ar9  computed  at  the  center 
of  each  square.  Wind-stress  and  bottom-friction  factors 
must  be  specified  or  computed  at  the  centers  of  each  square. 
This  staggered-grid  is  basic  to  the  spatial  realization  of 
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Figure  2.1  Staggered  Grid  o£  the  lumerical  Model 


Leendertse*  s  finite- differencing  schema  in  that  the  lean 
velocity  into  or  out  of  each  square  is  used  to  compute  the 
change  in  sea  level  within.  The  grid  also  permits  the 
coastal  constraint  of  zero  transport  perpendicular  to 
sea/land  boundaries. 

During  the  first  half  time  step,  implicit  computa¬ 
tions  proceed  row  by  row  from  left  to  right  and  explicit 
computations  proceed  column  by  column  from  the  bottom  to  th 
top  of  the  grid.  During  the  second  half  time  step  this 
process  is  reversed,  "centering"  tha  results  in  space  as 
well  as  time. 


2.  £21£&££  2I2SEL1  . 

The  computer  program  used  to  eodel  Honterey  Bay  is  a 
■odified  version  of  a  Fortran  prograe  developed  by  Hart 
(1976).  It  includes  provisions  for  nodeling  wind  stress, 
steady  flov  at  boundary  points,  and  overflow  at  boundary 
points,  none  of  which  were  available  in  Leendertse's  orig¬ 
inal  listing.  To  increase  the  flexibility  of  the  program, 
additional  modifications  have  been  made  during  this  study. 
These  include: 

•  Introduction  of  date  and  tine  computations 
to  permit  the  program  to  search  time-coded 
files  for  data  items  required  at  each  half 
tine  step. 

•  Direct  computation  of  bottom- friction 
factors  from  an  average  bottom- type  parameter 
or  from  a  grid  of  bottom  types. 

•  creation  of  an  interactive  subroutine  to 
start  the  program  by  prompting  the  user  for 
parameters  critical  to  each  run. 

•  Output  of  time  series  of  currents  and 
sea-surface  elevation  for  up  to  nine  points 
in  the  model  grid. 

•  Further  modularization  of  the  model's 
functional  components. 

The  core  of  the  program  is  subroutine  NODEL,  which 
contains  the  finite-differencing  scheme  of  Leendertse. 

Other  subroutines  serve  auxiliary  functions:  Interactively 
starting  the  run,  acquiring  both  constant  and  time-dependent 
data,  specifying  conditions  at  boundary  points,  specifying 
numerical  models  for  wind  and  bottom  stress,  and  supplying 
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results  in  various  output  formats.  An  outline  of  the 
prograa  is  presented  in  Table  I. 


TABLE  x 

Coaponents  of  the  Ruaerical  Model 

H AKEGRID 

Generate  depth  and  coaputation-control  grids. 

BATMODEL 

START 

DIMERS 

Main  prooraa  to  control  each  run  of  the  model. 
Interactively  input  rua-tiae  parameters. 

Input  constant  data  for  grid  of  area. 

FIND 

IRVAL 

INCOSB 

Locate  water  sections  to  be  modeled. 

Initialize  variables  everywhere  in  grid. 
Initialize  currents  as  desired. 

MODEL 

CHEZI 

BIND 

RESCLT 

Multioperational  f4.n4.te- difference  scheme. 
Supply  Chezy  coefficient  at  each  grid  point. 
Supply  wind  stress  tera  at  each  grid  point. 
Compute  output  values  in  desired  units. 

OPER 

STEAD! 

OVPLO 

OVFLD 

Specify  sea-surface  elevations  at  open  bounds. 
Specify  currents  at  open  bounds  and  rivers. 
Specify  overflow  currents  at  boundaries. 
Specify  overflow  threshholds. 

IRTIDE 
IN BIRD 

Get  tide  from  a  time-coded  file. 

Set  wind  from  a  time-coded  file. 

BEADS 

PRINT 

PLOT 

SERIES 

Output  headers  for  each  output  file. 

Output  results  for  desired  times. 

Output  results  needed  for  graphic  display. 
Output  results  for  specific  points. 

PTGRID 

CALERD 

ADTIHE 

Utility  print  of  input  gridded  data. 

Utility  for  number  of  days  in  month. 

Utility  to  increment  time. 

BAIPLOT 

Provide  graphic  presentation  of  results 

ELEVCCMP 

Compare  elevation  series  with  observed  data. 

COBRCOMP 

Compare  current  series  with  observed  data. 

Soae  subroutines,  such  as  OPEN,  aust  be  prepared 
specifically  for  the  area  to  which  the  nodel  is  applied, 
while  others,  such  as  MODEL,  should  not  be  altered. 
Hodularizatioa  of  the  prograa  peraits  the  user  to  readily 
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change  the  appropriate  functions  when  adapting  the  model  to 
different  coastal  areas  (and  computers)  . 

C.  USE  AND  ADAPTATION 

when  applying  the  numerical  modal  to  a  specific  area, 
certain  requirements  concerning  input  data  and  program  modi¬ 
fication  aust  be  set.  These  are  discussed  below. 

1-  l££fit  2&ta  £eflai£iaS£ts 

Each  cun  of  the  nuaerical  model  requires  specifica¬ 
tion  of  start  and  end  tines,  time-step  length,  and  an 
interval  at  which  results  aust  be  output.  ?or  experimental 
(as  opposed  to  operational)  use,  other  things  may  be  speci¬ 
fied:  Points  at  which  series  output  is  desired,  the  type  of 
output,  and  oaission  of  certain  terns  in  the  hydrodynamic 
equations. 

Input  values  that  are  unique  to  the  area  to  be 
aodeled  and  that  do  not  change  froa  one  run  to  the  next  must 
also  be  supplied.  These  values  ara  most  conveniently  stored 
in  a  separate  file.  They  include: 

•  A  location  title  and  central  latitude. 

•  Dimensions  of  the  grid. 

•  Depths  for  each  grid  corner. 

•  Control  numbers  for  each  grid  square 
(land*0,  water*  1,  overflow*2). 

•  A  general  bottoa-f r iction  parameter  or  a 
bottoa-type  indicator  for  the  center  of  each 
grid  square. 

•  Number  of  tide  stations  supplying  data  for 
boundary  conditions. 

•  Number  of  wind  stations  supplying  data. 
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Selection  of  the  size  of  grid  to  be  used  is  limited 
by  the  size  of  the  area  to  be  modeled  and  the  available 
virtual  memory  and  CPU  time  of  the  computer.  Since  depths, 
friction  and  find  factors,  output  data,  and  two  half-time- 
step  values  for  each  velocity  component  and  the  sea  surface 
elevation  must  be  available  for  each  grid  square  at  all 
times,  at  least  12  arrays  must  be  dimensioned  according  to 
the  grid  size  and  survey  area.  The  maximum  dimensions  of  80 
by  80  used  in  this  study  required  close  to  1  megabyte  of 
virtual  computer  memcry  and  approximately  0.5  s  of  CPU  time 
per  time  step  on  the  IBM  3033  mainframe  computer. 

2.  Subroutine  Modifications 

Stresses  at  the  bottom  and  surface  must  be  modeled 
in  the  subroutines  CHEZ?  and  WIND.  Since  Leendertse»s  model 
already  assumes  a  quadratic  formulation  for  bottom  friction, 
only  the  Chezy  factor,  C,  need  be  provided  by  the  subroutine 
CHEZY;  however,  many  empirical  techniques  exist  for 
computing  the  factor,  most  of  which  rely  on  a  description  of 
the  bottom  type.  The  user  may  select  the  Technique  which 
best  applies. 

Similarly,  the  user  must  program  a  wind-stress 
formulation  in  the  subroutine  WIND.  Values  for  wind  speed 
and  direction  are  obtained  from  time-coded  data  sets  using 
the  subroutine  INWIND. 

Subroutines  OPEN,  STEADY,  07 FLO ,  and  OVFLD  supply 
time-varying  values  for  sea-surface  elevation,  currents,  and 
overflow  conditions  at  both  open  and  closed  boundary  points 
in  the  grid.  In  OPEN,  an  algorithm  must  be  provided  to 
compute  the  variation  of  sea-surface  elevation  along  the 
open  boundaries  of  the  grid.  The  necessary  tidal  amplitudes 
are  obtained  from  time-coded  sets  of  data  for  established 
tide  stations,  using  subroutine  INTIDE.  STEADY  permits 
currents  to  be  assigned  to  individual  grid  points. 


overriding  the  computed  currents,  in  initial  current  field 
may  be  entered  using  subroutine  INCCJRfi.  Finally,  OVFLO  and 
OVFLD  permit  conditions  of  flooding  to  be  ascertained  and 
modeled  at  grid  points  assigned  the  control  value  2. 

ffhen  implementing  the  nodel  on  various  computers 
and  for  various  purposes,  modifications  may  be  necessary  in 
the  input/output  subroutines  START,  DINENS,  HEADS,  PRINT, 
PLOT,  and  SERIES.  RESULT  may  also  be  modified  to  compute 
additional  quantities  of  interest,  suca  as  horizontal 
transports. 

3 .  Computer  Implementation 

For  this  study,  the  numerical  model  was  implemented 
on  an  IBM  3033  mainframe  computer  at  the  Naval  Postgraduate 
School.  Only  a  few  minor  language  changes  were  required 
before  the  system's  Fortran  H  compiler  could  be  used  on  the 
program  originally  supplied  by  Hart.  Subsequent  modifica¬ 
tions  were  made  and  all  jobs  were  run  from  remote  terminals 
under  the  School's  interactive  time-sharing  system. 

The  system  made  possible  the  development  of  several 
programs  that  facilitated  preparation  of  data  for  input  to 
the  model  and  production  of  graphic  output.  Especially 
useful  among  these  were:  H AKEGRID,  a  program  that  generates 
input  depth  and  computation-control  grids  from  digital 
bathymetric  data  already  available  for  the  area,  simplifying 
the  otherwise  laborious  task  of  creating  a  grid  on  chart 
overlays;  programs  that  generate  predicted  tidal  amplitudes 
from  constituents,  or  from  data  supplied  in  the  NOS  Tide 
Tables;  and,  ELEVCOHP  and  CURRCOHP,  programs  that  plot 
time-series  output  from  the  model  against  observed  data  from 
the  same  time  period  for  verification  of  model  accuracy. 
Although  the  programs  themselves  may  not  be  transferable  to 
other  computers,  supplying  similar  auxiliary  software 
together  with  the  model  (or  even  incorporating  the 


algorithas  into  the  model  program)  greatly  enhances  the 
ready  application  of  the  aodel  to  other  coastal  areas. 

Data  necessary  to  running  the  nuaerical  model  were 
stored  in  computer  files  distinguished  by  type.  All 
constant,  gridded  data  were  stored  in  one  file  while  rime- 
varying  data  were  stored  in  separate  files  by  type  and  year 
(for  example,  MONTEREY  TIDE76).  In  this  way,  a  new  file  cf 
input  data  did  not  have  to  be  creatad  for  each  run  cf  the 
model.  The  sources  and  selection  of  input  data  are 
discussed  in  sections  3.  B  and  3.C. 


III.  APPLICATION  OF  JHE  MODEL  TO  MONJEHE Y  JAY 


A.  VALIDITY  OF  ASSBHPTIOMS 

The  applicability  of  the  numerical  model  to  Monterey  Bay 
was  checked  by  examining  the  assumptions  outlined  in  section 
2. A.  The  assumption  of  negligible  vertical  velocity  was 
confirmed  for  tidal  forcing  by  noting  that  the  maximum  depth 
of  the  area  modeled  is  much  less  than  the  wavelength  of  the 
semidiurnal  tide  (3  km  <<  7  600  km).  Other ,  horizontal, 
forcing  conditions  of  currents  and  wind  were  applied  as 
steady-state  phenomena  in  the  model.  In  addition,  since 
this  study  concerns  large-scale  fluid  motions  over  time 
periods  of  15  minutes  or  more,  short-period  turbulent 
effects  and  vertical  accelerations  were  neglected:  the 
hydrostatic  approximation  holds. 

Although  Monterey  Bay  is  not  vertically  well-mixed 
[Scott,  1973],  in  depths  of  a  few  hundred  meters  or  less  the 
difference  in  dynamic  height  between  that  of  the  assumed, 
homogeneous  density  profile  and  that  of  a  more  typical 

profile  is  less  than  1  cm,  which  is  negligible  for  the 
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purpose  of  this  study.  In  depths  of  1000  m  or  more,  the 
effect  is  more  significant  (several  centimeters)  ;  however, 
since  such  depths  occur  outside  the  bay  proper,  the  effect 
of  density  stratification  was  ignored  and  horizontal  veloci¬ 
ties  were  averaged  over  depth  to  obtain  a  general  picture  of 
circulation  in  the  bay. 

B.  COSSTANT  INPOT 

The  numerical  model  was  applied  to  three  different  grids 
covering  Monterey  Bay  (Figure  1.1  and  Figures  3.1  -  3.3). 
Grid  A,  a  small-scale,  1-  or  2-km  grid,  80  by  80  km,  was 


designed  to  perait  the  introduction  of  offshore,  non-tidal 
currents  as  a  steady  forcing  condition  and  to  place  the  bay 
far  enough  away  from  the  three  open  boundaries  to  remove 
their  associated  spurious  effects.  Since  this  grid  was 
particularly  vulnerable  to  nuaerical  instabilities  in  the 
model,  two  other  grids  were  devised.  Grid  B,  a  large-scale, 
1-lca  grid,  23  by  50  km,  covered  the  bay  and  reduced  the 
number  of  open  boundaries  to  one.  Grid  c,  a  1-km  grid,  40 
by  72  km,  covered  both  the  bay  and  sufficient  area  to  permit 
offshore,  non-tidal  currents  to  be  introduced.  All  grids 
were  skewed  20°  east  cf  north  to  align  the  boundaries 
perpendicular  to  the  tidal  forcing  conditions.  The  dimen¬ 
sional  and  constant  data  incorporated  into  these  grids  are 
discussed  below. 

1 •  Time  and  space  Dime nsions 

A  time  step  of  one  hour  was  chosen  to  permit  assess¬ 
ment  of  the  model  over  periods  of  several  days  without 
necessitating  extensive  use  of  CPD  time,  normally  the  model 
should  run  for  12-24  hours  (one  tidal  cycle  or  more 
depending  on  the  tidal  phase  differences  between  various 
parts  of  the  areal  to  establish  realistic  conditions  of 
current  and  sea-surface  elevation  throughout  the  area  [Hart, 
1976  ].  The  one-hour  interval  provided  a  sufficient  number 
of  data  values  fcr  comparison  with  hourly  or  half-hourly 
observations  of  sea- surface  elevation  and  currents. 

Since  computed  values  are  offset  in  the  staggered- 
grid  scheme  of  the  model,  use  of  a  relatively  small  grid 
size  in  regions  of  steep  bathymetric  relief  is  important  to 
model  accuracy.  With  the  constraint  on  array  dimensions  in 
mind,  the  smallest  grid  size  possible  (1  km2)  was  generally 
selected  to  permit  the  greatest  spatial  resolution  for  the 
area  of  concern. 
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Figure  3.1  Depth  contours  for  Grid  k 
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Figure  3.2  Depth  Contours  for  Grid  B 


32 


kuuiaiiWii 


X  (KM) 


Figure  3.3  Depth  coutoucs  for  Grid  c. 
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2.  jalkast£z  aaj  sain® 


Digital  bathymetry  was  provided  by  the  NOAA  National 
Geophysical  Data  Center,  where  depth  data  from  past  hydro- 
graphic  surveys  are  archived  for  most  coastal  areas  of  the 
United  states.  The  depths  were  positioned  by  latitude  and 
longitude  in  a  36-sec  grid  and  were  referenced  to  a  mean- 


lower-low-watar  (ML Li)  datum. 
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Figure  3.4  Locations  of  Bathymetric  Data  in  Grid  A. 


To  computer-generat  e  a  grid  of  depths  for  the  model, 
the  program  HAKEGRID  was  prepared.  The  program  first 
projects  the  bathymetric  data  onto  the  grid  coordinate 
system  (for  example,  grid  A),  which  is  a  modified-transverse 
Mercator  projection  skewed  about  a  specified  origin 


(Figure  3.4).  Depths  are  then  interpolated  at  the  corners 
of  each  grid  square. 

Since  depths  are  referenced  to  MLLW  and  a  straight¬ 
forward  correction  to  mean  sea  level  is  not  possible,  the 
datum  for  sea-surface  elevations  coaputed  by  the  model  was 
taken  to  be  MLLH.  All  input  tidal  amplitudes  were  likewise 
referenced  to  HLLH. 

Froa  the  gridded  depths  a  computation-control  grid 
wa3  automatically  generated  by  assigning  I's  to  all  water 
sguares  and  0's  to  all  land  squares  (assigned  dummy  eleva¬ 
tions  in  the  depth  grid)  .  Both  grids  could  be  altered  if 
necessary  before  their  use  in  the  model. 

3 .  Bottom  Friction 

To  model  bottom  stress,  the  empirical  Manning  equa¬ 
tion  for  the  Chezy  coefficient,  C,  was  used: 

C  »  (h+n)1/6/M  (3.1) 

The  coefficient  is  a  function  of  depth,  h,  instantaneous 
sea-surface  elevation,  n  ,  and  the  Manning  factor,  M,  which 
describes  bottom  roughness.  M  increases  with  bottom  rough¬ 
ness.  Clean  and  straight  natural  river  channels  typically 
require  H»0.025  to  0.030  m/s,  while  winding  channels  may 
require  H*0.033  to  as  high  a  value  as  0.15  m/s  in  very 
weedy,  overgrown  areas  [  p.  99,  Henderson,  1966]. 

Although  bottom  stresses  may  be  modeled  as  a  func¬ 
tion  of  the  bottom  type  or  texture  in  each  grid  square  (thus 
requiring  input  of  a  bottom-type  grid),  this  option  was  not 
exercised  for  Monterey  Bay.  over  the  large  area  covered, 
the  variation  in  depth  from  square  to  square  is  likely  to 
have  more  influence  than  the  relatively  small  variations  in 
bottom  roughness  that  occur  in  Monterey  Bay.  Following 
Spaulding  and  Beauchamp's  study  of  a  coastal  sea  (1983)  and 
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after  soae  exoer iaentation  (sectioa  IV. k),  a  constant  value 
of  H*0.04  a/s  vas  used  throughout  ths  bay. 

C.  TIHE-NABYING  INPOT 

The  forcing  conditions  of  tide  and  wind  were  accessed 
froa  time-coded  files.  Tidal  aaplitudes  were  applied  only 
along  the  open  boundaries  of  the  model,  whereas  wind  stress 
was  applied  over  the  entire  grid.  The  sources  and  applica¬ 
tion  of  these  data  are  discussed  in  this  section. 

1 .  Initial  &&d  itlojis 

At  the  start  of  a  run  of  the  model,  sea-surface 
elevations  were  approziaated  by  assigning  a  coastal  tidal 
amplitude  at  the  starting  time  to  every  point  in  the  grid. 
The  tidal  amplitude  at  Honterey  was  used  for  this  purpose. 
This  approximation  is  suitable  for  Monterey  Bay  since  the 
narrow  continental  3helf  and  the  absence  of  any  barrier 
islands  permit  the  tidal  wave  to  propagate  relatively 
rapidly  through,  the  area. 

A  zero  velocity  was  initially  assigned  to  each  grid 
point,  except  for  runs  including  an  offshore,  not-tidal 
current;  in  these  cases,  the  steady-state  current  velocity 
was  initially  assigned  to  offshore  grid  points. 

2.  Boundary  Hlal  Amplitudes 

A  major  factor  in  the  successful  application  of  the 
nunerical  model  vas  the  provision  of  suitable  tidal  forcing 
conditions  along  the  open  boundaries.  Tidal  aaplitudes  are 
predicted  in  the  NOS  Ii&e  Tables  J976  for  four  stations  near 
and  within  Honterey  Bay;  Ano  Nuevo,  Santa  Cruz,  Monterey, 
and  Carmel.  Tidal-constituent  aaplitudes  and  phases  are 
available  for  Honterey  and  Boss  Landing  (Appendix  A)  .  Since 
1963,  continuous  observations  of  water  level  have  been  made 
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at  Monterey.  A  two-year  series  of  nearly  continuous  obser¬ 
vations  was  aade  at  Moss  Landing  from  1976  through  1977.  To 
obtain  aaplitudes  along  the  boundaries  of  the  model  grid, 
coastal  values  such  as  these  must  be  extrapolated. 

Several  attenpts  have  been  aade  to  formulate  the 
effects  of  continental  shelves  on  the  open-ocean  tide 
[Clarke  and  Battisti,  1980;  Gill  and  Porter,  1980;  Munk,  et 
al. ,  1970].  Bacause  of  the  narrow  continental  shelf  and 
bisecting  canyon,  Monterey  Bay  does  not  satisfy  the  assump¬ 
tions  necessary  to  apply  these  foraulations.  However,  to 
gain  sons  insight  into  the  effect  of  the  extreee  depth 
difference  between  the  tide  station  at  Monterey  and  the 
offshore  boundary  points  of  the  model,  a  comparison  was  made 
between  tidal  constituents  obtained  at  Monterey  and  at  a 
pressure  gage  located  in  3903  m  of  water  offshore 
[Cartwright,  et  al. ,  1979].  The  results  are  presented  in 
Table  II. 

The  coastal  and  pelagic  phases  clearly  do  not  corre¬ 
spond  since  the  pelagic  gage  was  located  at  seme  distance 
from  Monterey  (see  Figure  1.1),  but  the  agreement  of  the 
amplitudes  suggests  that  the  aforementioned  depth  difference 
has  little  effect.  In  the  absence  of  any  more  certain 
method  for  extrapolating  tidal  amplitudes,  the  values  at  the 
coastal  station  were  applied  directly  along  a  line  of 
constant  phase  extending  out  from  shore. 

Examination  of  cotidal/cophase  charts  by  Munk,  et 
al.  (1970),  Luther  and  Wunsch  (1975),  and  Parke  and 
Henderschott  (1980)  ,  suggests  that,  in  the  vicinity  of 
Monterey  Bay,  the  tidal  wave  propagates  nearly  parallel  to 
the  coast.  The  model  grids  were,  therefore,  skewed  in  such 
a  way  that  the  open  boundaries  were  perpendicular  or 
parallel  to  the  coast.  At  each  time  step  the  forcing  ampli¬ 
tude  was  made  to  vary  directly  with  the  tide  at  Monterey  all 
along  the  southern  boundary,  with  the  tide  at  Ano  Nuevo  all 
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along  the  northern  boundary,  and  with  linearly  interpolated 
values  between  the  two  along  the  western  boundary. 

The  tidal  values  used  to  sat  the  boundary  conditions 
■ay  be  interpolated  from  the  NOS  Tide  Tables  1976,  computed 
from  constituents  (  Schureman,  1940],  or  taken  directly  from 
observed  data.  The  last  was  preferrl  since  the  first  two 
predictive  techniques  cannot  take  intr-  account  atmospheri¬ 
cally  forced  or  anoaalous  changes  in  se**  level,  such  as 
storm  surge.  However,  as  mentioned  previously,  observations 
were  available  only  at  Monterey  and  doss  Landing.  Some 
experimentation  was  necessary  to  model  the  phase  lags 
between  Monterey,  Santa  Cruz,  and  Ano  Nuevo  (section  IV. A). 
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3.  Boundary  Currents 


Two  types  of  flow  nay  be  imposed  at  the  boundaries 
of  the  numerical  model.  Pirst,  the  discharge  of  rivers 
along  otherwise  closed  boundaries  may  be  represented  as  a 
vertically  averaged  current  velocity  assigned  to  the  appro¬ 
priate  coastal  grid  point  during  each  time  step.  The  mean 
annual  discharge  of  all  major  streams  and  rivers  entering 
Monterey  Bay  is  1.85  x  10*  m3/day  [Broenkow  and  Smethie, 

1978  ],  which  amounts  to  a  vertically  averaged  current 
velocity  of  only  0.2  cm/s  were  all  rivers  to  enter  at  one 
point.  As  a  result,  the  river  inflow  was  judged  negligible 
for  this  application. 

A  second  type  of  flow,  currents  due  to  non-tidal 
effects,  may  be  imposed  in  offshore  regions  of  the  model. 

The  narrowness  of  the  continental  shelf  near  Monterey  Bay 
leaves  the  bay  particularly  open  to  forcing  by  large-scale 
oceanic  currents.  Previous  studies  of  the  area  suggest  that 
such  currents  are  an  important  force  driving  the  circulation 
of  the  bay  [Lazanoff,  1971;  Garcia,  1971;  Bretschneider , 
1982].  The  presence  of  an  offshore  current  was  simulated  by 
assigning  initial  velocities  to  the  offshore  portion  of 
grids  A  and  c  for  some  runs  of  the  model.  The  convective- 
inertia  terms  in  the  numerical  model  propagate  the  current 
influence  into  the  inshore  portions  of  the  grid.  A  north¬ 
ward  current  of  25  cm/s  (0.5  knots)  was  assigned.  This 
value  was  proposed  by  Scott  (1973)  as  a  simple,  steady-state 
aodel  for  the  offshore  circulation. 

4.  Wind 

Wind  stress,  F^,  was  parameterized  within  the  numer¬ 
ical  model  by  the  familiar  quadratic  law  [Dronkers,.  1964]: 


F. 

1 


(3.2) 


=  <2.6x10~3)p  w|w|/p(h+n) 

SL 

W  represents  the  wind  velocity  vector,  p,  is  atmospheric 
pressure,  p  is  the  mean  density  of  the  water,  h  is  the 
depth,  and  n  is  the  time-varying  sea-surface  elevation.  The 
model  permits  input  of  wind  direction  and  amplitude  as  a 
forcing  condition  over  the  whole  field  of  the  grid  for  a 
specified  range  of  time  steps.  Monthly  distributions  of 
wind  at  Santa  Cruz  and  Moss  Landing  for  the  period  May, 

1976,  through  May,  1977,  were  obtained  from  the  Santa  Cruz 
Wastewater  Facilities  Planning  Study  [Brown  and  Caldwell, 
Inc.,  1978],  Average  and  maximum  values  for  the  wind  were 
applied  to  some  runs  of  the  model  (see  section  IV. B)  . 

0.  DATA  FOB  COMPARISON 

The  numerical  model  was  calibrated  by  comparing  modeled 
sea-surface  elevations  and  current  velocities  at  specific 
grid  points  with  observed  values  at  the  same  locations.  The 
comparison  process  was  limited  by  a  paucity  of  suitable, 
long-term  water-level  and  current-meter  observations  for 
Monterey  Bay.  Water-level  data  are  available  only  fcr 
National  Ocean  Service  tide  stations  at  Monterey  and  Moss 
Landing.  The  primary  sources  for  current-meter  data  are 
predesign  studies  conducted  for  the  emplacement  of  sewage 
outfalls  near  Santa  Cruz,  the  Pajaro  River,  and  the  Salinas 
River,  but  only  data  for  Santa  Cruz  and  the  Salinas  River 
could  be  obtained.  In  some  fortuitous  instances,  both 
water-level  and  current-meter  data  were  collected  concur¬ 
rently  (Figure  3.5)  .  The  period  July-August  1976  yielded 
sufficient  data  for  comparisons  at  the  two  water-level  and 
at  three  current-meter  stations;  their  locations  are  plotted 
on  each  model  grid  (Figures  3.1  -  3.3). 
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Figure  3.5  Observing  Periods  for  Comparison  Data. 

1  •  Water-level  Observations 

Water-level  observations  have  been  made  nearly 
continuously  since  1963  at  NOS  Tide  Station  9  41-3450  on  the 
seaward  end  of  Municipal  Wharf  2  in  Monterey.  The  float- 
type  tide  gage  is  located  in  8.2  m  of  water.  Recorded  times 
are  accurate  to  within  6  minutes  and  heights  are  resolved  to 
3.0  cm  [  Bretschn eider,  1980  ].  Digitized  hourly  heights  for 
the  period  11/8/73  to  3/2/83  were  obtained  from  the  National 
ocean  Service,  Tidal  Datums  Section  N/OMS123.  The  heights 
were  corrected  to  MLLW  and  converted  from  feet  to  meters. 

In  addition  to  providing  comparison  data,  these  observations 
were  used  to  determine  boundary  amplitudes  for  some  runs  of 
the  model. 

At  Moss  Landing  water-level  observations  were  made 
for  20  months  as  part  of  the  California  Marine  Boundary 
Program  [National  Ocean  Survey,  1981].  NOS  Tide  Station 
94  1-3616  was  a  float-type  tide  gage  located  at  the  seaward 
end  of  the  Moss  Landing  Ocean  Pier  in  9.1  m  of  water. 
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Digitized  hourly  heights  were  obtained  for  the  entire  period 
5/9/76  to  1/10/78  and  processed  as  for  the  station  at 
Monterey. 

2 .  Current- m^gr  Obser vations 

As  part  of  the  Santa  Cruz  Wastewaters  Facilities 
Planning  Study,  a  current-o eter  station  was  located  1  mile 
offshore  of  Terrace  Point  in  30  a  of  water  [Brown  and 
Caldwell,  Inc.,  1978  ].  Two  AMF  Vector  Averaging  meters  were 
installed  at  9-  and  15-m  depths  for  the  periods  June  to 
November,  1976,  and  January  to  May,  1977.  The  only  data 
that  could  be  obtained  for  comparison  purposes  covered  July 
and  August,  1976,  at  the  15-m  depth.  The  data  included 
7.5-ainute  averages  of  current  speed  and  direction, 
expressed  as  a  pair  of  orthogonal  velocity  vectors. 

Two  current-meter  stations  were  occupied  approxi¬ 
mately  1  nautical  mile  north  and  south  of  the  Salinas  River 
during  oceanographic  investigations  for  the  Monterey 
Peninsula  Water  Pollution  Control  Agency 

[Engineering-Science,  Inc.,  1977  ].  At  each  station,  two 
ducted-impeller- type  meters  were  installed  at  9  and  15  m  for 
the  overall  period  January,  1976,  to  January,  1977.  Current 
speeds  and  direction  were  averaged  at  30-minute  intervals 
and  expressed  as  a  pair  of  orthogonal  velocity  vectors. 
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A.  COMPARISON  OP  BO DEL  RESULTS  WITS  OBSERVATIONS 

Comparison  of  the  model  with  observations  was  used  both 
to  "fine  tune,"  or  calibrate,  the  numerical  model,  and  to 
assess  its  general  validity.  The  affects  of  varying  input 
constants  such  as  the  size  and  resolution  of  the  grid,  the 
time  step,  and  the  bottom-f riction  coefficients  were  consid¬ 
ered.  In  addition,  schemes  for  detarmining  boundary  tidal 
amplitudes  and  for  including  non-tidal  current  fields  were 
tested  in  an  affort  to  match  observed  elevations  and 
currents  as  closely  as  possible. 

Application  of  the  model  to  grid  A  revealed  apparent 
numerical  instabilities  that  caused  overflow  in  the  computa¬ 
tions  after  as  few  as  1.25  days  (33  time  steps).  The  sudden 
oscillations  in  sea-surface  elevation  at  Honterey  were  due 
to  propagation  into  the  bay  of  extreme  amplitudes  and 
currents  generated  in  the  offshore  portion  of  the  grid 
(Figure  4.1).  The  overflow  condition  was  unaffected  by 
changing  the  resolution  of  the  grid  from  2  to  1  km,  but  was 
very  sensitive  to  changes  in  the  phasing  of  the  tidal  ampli¬ 
tudes  along  the  open  boundaries.  Under  the  premise  that  the 
presence  of  three  open  boundaries  enhanced  instabilities, 
grid  B  (with  one  open  boundary)  and  grid  C  (with  two  open 
boundaries)  were  subsequently  used  during  the  comparison 
process. 

1-  HSZaUoa  *02£a£isoi}s 

A  run  of  the  model  for  twelve  days  at  the  one-hour 
time  step  produced  seme  agreement  between  modeled  and 
observed  sea-surface  elevations  at  .lonterey  and  very  good 
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Figure  4.1  sea-surface  Elevations  at  Monterey,  Grid  A. 

agreement  at  loss  Landing  (Figures  4.2  and  4.3).  Varying 
the  Manning  bo  tt  om- friction  factor  improved  this  result 
(Table  3.1  and  Figures  4.4  and  4.5). 


TABLE 

Ill 

Effect  of  Various 

Banning 

Factors 

M  (m/s)  =  .03  .04 

.05 

.06  .10 

Monterey  7.4  6.8 

6.3 

5.8  4.8 

Moss  Landing  4.6  4.5 

4.4 

4.3  4.0 

Table  values  are  the  RMS  errors  in 

centimeters 

for  comparisons  between  modeled  and 

observed 

sea-surface  elevations. 

j 

The  seemingly  unrealistic  Manning  factor  of  0.10  m/s 
■ay  serve  as  a  description  of  the  bottom  roughness  over  the 
large  area  of  each  grid  square,  if  roughness  is  thought  of 
in  terms  of  the  steep  slope  that  is  otherwise  not  accounted 
for  in  the  model.  The  higher  Manning  factors  did  not, 
however,  improve  the  results  of  current  comparisons;  they 
served  primarily  to  damp  out  noise. 


Figure  4.2  Elevations  at  Honterey,  At=3600  s. 


An  oscillation  that  appeared  forced  by  the  time  step 
was  evident  in  the  modeled  curves.  It  was  especially 
evident  at  HOQterey  and  when  a  shorter,  15-minute  time  step 
was  used  (Figures  4.6  and  4.7). 

The  noise  may  be  the  result  of  applying  observed 
tidal  amplitudes  as  the  boundary  forcing  condition.  The 
observed  water  levels,  digitized  hourly,  may  include  jumps 
and/or  contaminating  frequencies  due  to  the  recording 


Figure  4.3  Elevations  at  Boss  Landing,  At=3600  s 


I 


instruments.  Linear  interpolations  required  between  the 
hourly  amplitudes  for  each  half  time  step  may  have  exagger¬ 
ated  the  instrumental  effects.  To  better  judge  rhe  use  of  a 
shorter  time  step,  raw  water-level  observations,  usually 
made  at  a  6-minute  interval,  should  be  applied  to  the  model. 

Spectral  analysis  of  the  modelad  and  observed 
curves,  in  addition  to  reflecting  their  general  agreement, 
reveals  spurious  frequencies  generated  by  the  model  at 
Monterey  (Figures  4.8  and  4.9).  The  incoherent  frequencies, 
which  are  also  found  in  the  spectra  for  currents  at  Santa 
Cruz  (Figure  4.14),  correspond  to  apparent  periods  of  3.0 
and  2.2  hours.  These  periods  are  longer  than  the  1-hour, 
fundamental  seiche  period  for  Monterey  Bay  [Lynch,  1970  ]. 
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Figure  4.5  Elevations  at  Boss  Landing,  M*. 10  a/s. 
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igure  4.6  Elevations  at  Bontarey,  At*900  s 
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Figure  4.7  Elevations  at  Boss  Landing,  A  t«900 
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Pigure  4.9  spactra,  Elevations  (SE)  at  Boss  Landing. 
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The  phenomenon  of  time-step-linked  oscillations  has 
been  experienced  by  other  investigators  [Chiang  and  Lee, 
1982],  who  found  that  generating  a  latheaatically  saooth 
function  froa  the  observed  data  provided  more  suitable 
amplitudes  for  forcing  the  model.  A  smooth  amplitude  func¬ 
tion  was  obtained  in  this  study  by  summing  the  tidal 
constituents  for  Monterey.  Applying  boundary  conditions 
based  on  these  predicted  tides  gradually  reduced  the  noise 
in  the  model  (Pigure  4.10).  Mot  surprisingly,  however,  the 
predicted  elevations  did  not  provide  good  agreement  with  the 
observed  water-levels. 


o 


Pigure  4.10  Elevations  Using  Predicted  Boundary  Amplitudes. 

The  sensitivity  of  the  model  to  open- boundary  condi¬ 
tions,  already  noted  in  the  case  of  grid  A,  was  demonstrated 
by  comparing  results  for  different  phase  lags  along  the  one 
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boundary  of  grid  B  (table  IV)  .  k  constant  amplitude  along 
tha  boundary  (no  phase  lag)  was  found  to  be  nost  successful, 
though  the  HOS  Tables  1976  predict  that  Honterey  lags 

Santa  Cruz  by  6  minutes. 


TABLE  IV 

Effect  of 

Boundary  Amplitude  Phasing  on  Srid  B 

Phasing 

6+  0  6- 

Honterey 

Boss  Landing 

7.9  7.4  7.8 

5.3  4.6  5.0 

Table  values  are  the  SMS  errors  in  centimeters. 
Phasing  is  Honterey  minus  Santa  Cruz,  in  minutes. 

2 .  Current  Coaparison3 

The  model-generated  current  values  agreed  poorly 
with  the  15-m  depth  observations  at  all  three  locations 
where  coaparisons  were  made  (Figures  4.11-4.13).  In  making 
the  comparisons,  current  vectors  from  the  aodel  output  and 
froa  the  observed  records  were  resolved  into  eastward  and 
northward  components,  taking  into  account  their  respective 
skewed  coordinate  systems.  The  modeled  current  components, 
particularly  near  the  Salinas  River  mouth,  were  weak  or 
non-existent.  Clearly,  forces  in  addition  to  tides  were  at 
work  in  generating  the  observed  currents. 

For  the  currents  near  Santa  Cruz  the  apparent  rough 
correspondence  in  frequency  was  confirmed  somewhat  by  a 
spectral  analysis  (Figure  4.14).  The  gradual  increase  in 
the  northward  component  of  the  modeled  current  may  reflect 
long-period  variations  in  the  current  field. 


In  an  effort  to  obtain  better  current  agreement  ard 
to  see  the  effects  of  currents  associated  with  the  conti 
neatal  slope,  a  steady  offshore  current  was  applied  to  grid 
C.  The  results  near  Santa  Cruz  were  poor;  inherent  oscilla¬ 
tions  in  the  lodel  were  amplified  and  agreement  was  not 
improved.  The  model's  sensitivity  to  open-boundary  condi¬ 
tions  made  this  a  difficult  subject  to  pursue  within  the 
scope  of  this  study. 

B.  MODELED  CXRCOLATIOH  OF  HOHTEBET  BAT 

Although  the  tidally  forced  numerical  model  did  not 
reproduce  observed  currents  at  the  comparison  locations,  it 
should  nevertheless  provide  an  estimate  of  the  barotropic 
tidal  circulation  of  Monterey  Bay.  A  general  view  of  circu¬ 
latory  patterns  nay  be  obtained  by  examining  the  modeled 
sea-surface  elevation  and  current  fields  for  a  24-hour 
period. 

1*  liiaiil  £2£ce£  Simulation 

A  tidally  forced  elevation-field  series  is  presented 
in  Appendix  B.  A  small  oscillatory  structure  in  the 
southern  bight  is  consistent  with  the  greater  amplitude  of 

noise  observed  at  Monterey.  In  other  respects,  the  eleva¬ 

tion  fields  generally  show  a  clear  progression  of  the  tidal 
wave  into  the  bay. 

A  series  of  current-field  plots,  including  the  volu¬ 
metric  transport  associated  with  each  vector,  is  presented 
in  Appendix  C.  A  total  volume  of  2.0x10»  appears  to  be 

pumped  across  the  boundary  during  the  diurnal  tidal  cycle. 

This  barotropic  result  may  be  contrasted  with  the  10’  a3 
pumped  by  internal  tidal  mixing  reported  by  Broenkow  and 
Smethie  (1978)  . 
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The  current  plots  show  generally  weak  (<5  cm/s) 
tidal  flow  into  and  out  of  the  bay,  corresponding  to  periods 
of  flood  and  ebb.  In  the  southern  parr  of  the  bay,  a  strong 
east-west  current  (up  to  30  cm/s)  appears  just  north  cf  the 
Monterey  peninsula.  This  jet  is  consistent  with  the  strong 
currents  experienced  by  divers  in  the  area.1  In  the  northern 
part  of  the  bay,  a  broader  current  (10  cm/s)  flows  along  the 
depth  contours. 

Of  especial  interest  is  a  current  pattern  that 
develops  between  Ano  Nuevo  and  Santa  Cruz  on  current  plots 
made  using  grid  A  (Figure  4.15).  The  gyre,  which  the  model 
predicts  to  have  speeds  ranging  from  2  to  10  cm/s,  is 
consistent  with  the  observations  of  Carter  and  Kazmierczak 
(1963)  who  noted  a  closed  circulation  in  the  area  with 
similar  speeds. 

2.  T id ally  Forced  Circulation  with  Hind 

When  an  average  wind  of  3  m/s  (7  mph)  from  the 
west -northwest  was  applied  over  the  entire  field  of  grid  B, 
the  tidally  forced  circulation  was  unchanged.  A  maximum 
wind  of  10  m/s  (30  mph)  from  the  west-northwest  also 
produced  essentially  unchanged  circulatory  patterns.  A 
12-day  series  of  modeled  elevations  at  Monterey  and  Moss 
Landing  under  the  same  maximum  wind  yielded  results  nearly 
identical  to  those  without  wind. 


'Per  conversation  with  Dr.  E.  C.  Haderlie,  Dept,  of 
Oceanography,  Naval  Postgraduate  Schoox,  9/27/83. 


Figure  4.15  Current  Field  for  Srid  k  at  76  0  7  0  2  1500 
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V.  COBCLOSIOBS 


A.  VALIDITY  DP  THE  MODEL 

The  numerical  model  has  been  shown  by  comparison  with 
observations  to  provide  reasonably  accurate  sea-surface 
elevations.  Although  the  modeled  current  does  not  accounr 
for  the  total  observed  current  at  comparison  stations,  it 
probably  accurately  reflects  the  contribution  due  to  the 
barotropic  tide.  Much  of  the  remaining  observed  current  may 
represent  response  to  diurnal  wind  stress,  offshore  non- 
tidal  currents,  and/or  forcing  due  to  internal  waves,  the 
last  requiring  a  more  complex,  three-dimensional  model  for 
further  investigation.  Application  of  observed  winds  on  a 
timestep-by-timestep  basis  might  be  a  fruitful  avenue  for 
further  investigation. 

That  numerical  instabilities  exist  in  the  model  has  been 
noted  by  various  authors  [Soe,  et  al.,  1978;  Benque,  et  al. , 
1982  ],  who  have  proposed  some  improvements  in  the  model's 
formulation.  The  flexibility  and  accuracy  of  the  model 
might  be  improved  by  further  investigation  of  their 
suggestions. 

B.  HTDBOGBAPHIC  SURVEY  APPLICATIONS 

A  numerical  model  such  as  that  implemented  by  this  study 
can  improve  the  process  of  correcting  depths  for  changes  in 
sea-surface  elevation  during  a  hydrographic  survey. 
Advantages  of  the  two-dimensional  model  over  simple  extrapo¬ 
lative  techniques  or  more  complex,  three-dimensional  models 
result  from  the  model's  relative  simplicity,  flexibility, 
ability  to  operate  in  a  real-time  data  collection  system, 
and  ability  to  compute  sea- surface  elevations  with  suffi¬ 
cient  accuracy. 


VS-..-S 


The  tasted  aodel  is  a  relatively  simple  FORTRAN  program 
to  implement  and  use,  particularly  with  the  interactive 
modifications  made  during  this  study.  It  could  be  further 
improved  in  this  respect  by  adding  an  interrupt/restart 
routine  to  permit  changes  of  constants,  such  as  time  step  or 
output  frequency,  during  the  course  of  a  single  computer 
run.  The  model  is  flexible,  since  it  can  be  readily  applied 
to  various  coastal  areas  by  means  of  the  grilling  software 
developed  during  this  study. 

To  implement  the  model  on  a  microprocessor  during  data 
collection  in  the  field,  requires  the  availability  of  suffi¬ 
cient  virtual  storage  capacity  and  CPU  time  to  permit  unim¬ 
peded  computations.  The  virtual  storage  required  depends 
upon  the  dimensions  cf  the  grid;  a  more  economical  use  of 
arrays  in  the  aodel  program  can  reduce  the  requirement.  In 
the  real-time  mode  of  operation,  computations  should  immedi¬ 
ately  follow  boundary-amplitude  updates  at  each  half  time 
step  to  make  efficient  use  of  CPU  time.  At  the  conclusion 
of  each  full  tima  step  the  resulting,  updated  sea-surface 
elevations  are  then  promptly  available. 

The  time  step  used  depends  upon  the  interval  at  which 
water-level  observations  are  available  from  one  or  more 
locations  suitable  for  establishing  boundary  amplitudes.  In 
the  real-time  mode,  presumably  such  data  could  be  teleme¬ 
tered  to  the  survey  vessel  at  the  standard  tide  gage 
frequency  of  6  minutes,  permitting  a  model  time  step  of  12 
minutes.  Since  updated  elevations  can  only  be  available  at 
the  conclusion  of  a  time  step  (Figure  5.  1)  ,  a  6-minute  lag 
exists  that  may  be  removed  only  by  post-survey  processing. 

Another  factor  in  real-time  operation  cf  the  model  is 
the  start-up  time  required.  The  model  should  be  calibrated 
to  establish  the  validity  of  friction  models  and  boundary- 
amplitude  algorithms,  preferably  by  comparing  the  output  for 
several  tidal  cycles  with  observed  tides  at  a  location  in 
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Figure  5.  1  Time  Step  Lag  During  Real-time  Model 


the  interior  of  the  survey  area.  If  historical  data  are  not 
available,  this  could  reguire  several  days  of  observation 
and  analysis  prior  to  the  survey. 

The  accuracy  of  the  model  in  computing  sea-surface 
elevations  from  tidal  forcing  alone  has  been  estimated 
during  this  study  as  4-8  cm  (t  RMS  error).  Survey  require¬ 
ments  are  3a  <  9. 14  ♦  0.005h  cm,  where  a  is  the  standard 
error  and  h  is  the  depth  in  centimeters  [Mobley,  1982].  The 
depths  at  the  tide  gages  in  this  study  were  about  8  m, 
permitting  a  3a  aqual  to  13.1  cm.  This  value  was  attained 
at  Moss  Landing  (3  RMS  error  *  12  cm)  and,  if  the  trouble¬ 
some  noise  could  be  removed  by  post-survey  processing  of 
water-level  observations  to  obtain  a  smooth  tidal  forcing 
function,  it  may  be  generally  attainable. 

The  model  itself  requires  further  development,  both  in 
the  application  to  Monterey  Bay  and  in  general.  Further 
testing  of  the  application  to  Monterey  Bay  should  include 
the  introduction  of  time-varying  wind  and  oceanic  currants, 
as  well  as  forcing  at  a  time  step  using  tidal  amplitudes 
observed  at  a  shorter  interval  than  the  1-hour  interval  used 
in  this  study.  The  model  should  also  be  applied  to  and 


quantitatively  tested  against  other  coastal  configurations 
and  bathyaetries  to  ensure  that  sufficient  accuracy  can  be 
consistently  achieved. 
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